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We describe interacting bosons at low temperature in spatially correlated random potentials. By 
a Bogoliubov expansion around the deformed mean-field condensate, the fundamental Hamiltonian 
for elementary excitations is derived, achieving an analytical formulation in the case of weak disor- 
der. From this, we calculate the sound velocity and true zero-temperature condensate depletion in 
correlated disorder and all dimensions. 
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The interplay between interaction, quantum statistics, 
and randomness defines one of the richest problems in 
condensed matter physics: the dirty boson problem ^[2]. 
Determining the ground state of a disordered Bose gas 
is already a formidable task [3]; all the more desirable 
is an analytical theory for the excitations of disordered 
Bose-Einstein condensates (BECs). In presence of a well- 
developed condensate, i.e. at low temperature and for 
weak disorder, the most economic description, due to 
Bogoliubov [4], treats quantum fluctuations around the 
best mean-field approximation to the condensate. One 
of the key quantities that such a theory should provide 
is the experimentally measurable, disorder-renormalized 
sound velocity, which characterizes the low-energy dis- 
persion and enters all thermodynamic properties. In this 
respect, existing theories for bosons in disorder are not 
entirely satisfactory. While simple approaches cannot de- 
termine a change in quasiparticle dispersion at all [5H2], 
more elaborate calculations find a positive correction to 
the sound velocity due to uncorrelated disorder [8MlO|. 
Still others, in different settings and with different meth- 
ods, report a negative correction [TTMl4| . Here, a unify- 
ing framework is missing. On a more conceptual level, the 
condensate deformation by disorder is often confounded 
with the condensate depletion, i.e. the fraction of particles 
not in the condensate at all. This quantum depletion is a 
crucial quantity whose smallness validates Bogoliubov's 
approach; it contains all particles with non-zero momen- 
tum [15] only if the condensate is homogeneous. If the 
condensate is deformed, it has non-zero momentum com- 
ponents on its own, already on the mean-field level [16] . 
Only these have been counted by previous approaches 
[5HTU]. For the true disorder-induced quantum depletion 
there are, to our knowledge, only few numerical results 
|17j . An analytical calculation of this pivotal quantity is 
lacking so far. 

With this Letter, we present a comprehensive Bogoli- 
ubov theory for inhomogeneous Bose-Einstein conden- 
sates that gives access to a wealth of relevant quantities, 
including the full excitation dispersion with the renor- 
malized speed of sound and the localization length of 
elementary excitations [TBI US]- Notably, our formula- 
tion encompasses the class of spatially correlated disorder 



that is currently under study with ultracold gases | 
By expanding the many-body Hamiltonian around the 
deformed mean-field solution, we derive the fundamental 
Bogoliubov Hamiltonian for excitations, verified to be or- 
thogonal to the ground state. For weak disorder, a fully 
analytical description is achieved. From this Hamilto- 
nian, we calculate corrections to the sound velocity and 
zero-temperature quantum depletion for correlated dis- 
order in all dimensions. 



STARTING LINE 

Interacting bosons are described by 

E = Jd d r&(r) h(r) + ^(r)^(r) #(r), (1) 

with the grand-canonical single-particle Hamiltonian 

h{r) = -^V 2 + V{r)-v (2) 

and field operators that obey [^(t), wT(r')] = S(r — r') 
[T5] . The global confining potential is assumed to be very 
smooth, ideally a very large box, and V{r) describes lo- 
cal spatial fluctuations. Repulsive interaction between 
bosons is accounted for by g = ATrh 2 a s /m > (in d = 3, 
with s-wave scattering length a s ), an excellent approxi- 
mation for cold and dilute gases, where the gas parameter 
(na;?) 1 / 2 is small. 

Below a critical temperature, the Bose gas forms a 
BEC, i.e. a large fraction of particles condense into the 
ground state of the single-particle density matrix. In the 
absence of interaction, this is the ground state of the po- 
tential V(r). Also interacting bosons condense, into a 
mode whose shape results from the competition between 
kinetic energy, confinement and interaction. Bogoliubov 
theory [4] takes advantage of this macroscopic occupation 
and splits the quantum field into a mean-field condensate 
and quantized fluctuations: ^(r) = $(r) + 6^(r). 
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CONDENSATE 

We first describe how a weak external potential de- 
forms the condensate. By definition, the ground state 
minimizes the energy ([I]) on the mean-field level and 
thus obeys the stationary Gross-Pitaevskii (GP) equation 
/i(r)$(r) + g$(r) 3 = [TS]. The condensate's kinetic en- 
ergy is minimized by choosing a fixed global phase, and 
we can take $(r) real. The imprint of a weak potential 
on the condensate amplitude can be computed perturba- 
tively by expanding $(r) = (r)+<f>( 2 ' (r)+. . . in 

powers of V around the homogeneous solution <&(°) = sfn 
[TT1 ITS] , In order to maintain a fixed average density 
n, also the chemical potential is adjusted at each order, 
H = gn + /iW + /i^ 2 ) + . . . . We insert these expansions 
into the GP equation, transform to momentum repre- 
sentation, and collect orders up to V 2 . The first-order 
imprint then is 



$ (1) = (1- 



2<?n 



(3) 



This linear-response deformation is proportional to the 
potential's matrix element V q = L~ d J d d re~ lq r V(r). 
In the denominator of eq. ([3]), the bare kinetic energy 
Eg = h 2 q 2 /2m equals the interaction energy gn when q 
equals the BEC healing length £ = h/y/2mgn. Thus, 
the condensate readily follows potential components with 
5$ « 1, but shows a strongly smoothed imprint when 
<?£ 3> 1 |16) . Pushing the expansion to second order, we 
find 



*( 2 ) = _J_ $ (D ^(i) (1 - M<£ ~ g" 
<J jyi/2 A/ i-P p 2gn + e° 



(4) 



Eqs. ([3]) and Q determine the disorder imprints also on 
derived quantities like the density rik — L~ d ®k-q&q- 



FLUCTUATIONS 

We expand the Hamiltonian ([I]) in powers of S^> and 
<5^t around the condensate. The linear term vanishes, 
because $(r) minimizes the energy functional. The rel- 
evant contribution is then the quadratic part that can 
be readily expressed in density-phase variables <5\l/ = 
<5n(r)/2$(r) +i$(r)6(p(r) [23]: 



H 
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2m 



2$(r) 
i> 2 (r)(V6^) 2 



[V 2 $(r)] 
4$ 3 (r) 

+ |„ 2 



Sn 2 



(5) 



In a homogeneous system with <i>(r) = y/n, it is advis- 
able to transform to Fourier space and Bogoliubov exci- 
tations 

Sfik 



7fc 



2ai 



ia k VnSip k 



(6) 



These bosonic excitations obey [%,T k ,] = Skh'- More- 
over, by choosing a k = (efc/efc) 1 ^ 2 the homogeneous 
Hamiltonian becomes diagonal, 



# (0) - E tffa 



(7) 



with efc = [e ( l(2gn + e^)] 1 ' 2 the clean Bogoliubov disper- 
sion relation [1]. Characteristically, low-energy excita- 
tions (kt; <C 1) have the linear dispersion = ficfc, with 
sound velocity c = -J ' gn/m. 

In presence of disorder, two things change. First of all, 
the reference point, from where the excitations originate, 
is shifted to the inhomogeneous ground state $(r). But 
still, eq. ^ is fit to define excitations with the proper 
commutation relations. The density-phase representa- 
tion ([6]) implies also that the fluctuation 



<5*(r) = J2\- U k( r )% ~ «fc(r)*7fe] 



decomposes over excitations with modes 
Uk(r) 
Vk(r) = 



1 

2L d / 2 
1 



2L d / 2 



$0) aky/n ' 
dkVn $(r) 
$(r) a k \/n 



(8) 

(9) 
(10) 



For the clean system with $(r) = s/n, these modes 
reduce to plane waves with the well-known amplitudes 
Uk = + a k) an d v k — ^(a^ 1 — a k ) such that 

u 2 ~ v 2 . = 1. In presence of disorder, the modes ^ 



and ( 10 1 are defined such that they still satisfy the bi- 



orthogonality 



d d r [ul(r)uk'(r)-vl(r: 



Vk'Vr) 



5. 



kk' 



(11) 



required for eigenmodes of the Bogoliubov Hamiltonian. 
Moreover, the excitations are also orthogonal to the de- 
formed ground state, because $(r) [ufe(r) — Ufe(r)] is a 
plane wave with zero average for all k 7^ [24l |25] . These 
are crucial properties for the low-energy excitations of the 
system to be well defined |26j . 

As a second difference to the homogeneous case, these 
excitations now live on a deformed background. Both 
differences can be accounted for by defining a single ef- 
fective potential Vkk' — ( y' w)kk' me diates scatter- 
ing between the different components of the Bogoliubov- 
Nambu pseudo spinor Tk = (7fc, 7!^)'- From ([5]), we thus 
arrive at the inhomogeneous Bogoliubov Hamiltonian 



H 



(12) 



fe.fe' 



with the structure H = i? (0) + At this point, 

the only approximation made is the neglect of third and 
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fourth order terms in the fluctuations. In contrast, 
is still exact in the disorder strength. 

A perturbative, but fully analytical description is ob- 
tained by expanding V = ( ™ £ j = V (1) + V (2) + . . . to 
lowest orders in the bare disorder with the help of eqs. 
|3]) and Q. The small parameter of this expansion is 
v = V/gn <C 1. The first-order scattering amplitudes 

W kk> = Wkk>Vk-k> and Y$ = y^V h -w are propor- 
tional to Vk-k' , as required by conservation of momen- 
tum. All information about interaction and condensate 
background is factorized into amplitude envelopes 



w 



(!) _ a fc a fc <£ 2 (l - 6 kk >) 



kk' 



Vkk' 



2 + £ 2 (fc'- k) 2 
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(13) 



Second-order scattering amplitudes are later only needed 
for k — kl because the disorder average restores transla- 



tion invariance: W}^ = — 



= E 



(2) 
W, 
p kp 



Vk-pV, 



p—k 



with 



3fc 2 



(2) afc V + 3(fc-p)' , „.. _ 
^ fe f 2gn [2 + {k- P ye? ( P '' 1 ' 

This concludes our derivation of the inhomogcncous Bo- 
goliubov Hamiltonian. From here, one can derive nu- 
merous physical quantities for any given potential V(r). 
For notational simplicity only, we assume in the follow- 
ing that V(r) describes disorder that is homogeneous and 
isotropic under the ensemble average, with V = and 



V q V- 



% q .V 2 (j d C d {qcr). 



(15) 



The dimensionless function Cd{q<?) characterizes the po- 
tential correlations persisting on the length scale a; the 
normalization is chosen such that in the thermodynamic 
limit / 0^C d (u) = 1. 

LOCALIZATION LENGTH 



The Hamiltonian (12) is a random operator, varying 



with each realization of the quenched disorder poten- 
tial. Therefore, the Bogoliubov excitations in d = 1 
are expected to be localized by the disorder [T8J IT9"] . 
And indeed, we can calculate their localization length 
as 'ioc = 72fc/2i> g from the backscattering rate 72ft and 
group velocity hv g = d^k- To lowest order in the small 
parameter v = V/gn <C 1, the backscattering rate de- 
rived by Fermi's Golden Rule from the Hamiltonian ( 12 ) 
reads 



72fc 



(16) 



where the density of states is p(ek) — [irhVg\ 1 . The re- 
sulting = \v 2 k 2 aCi{2ka)/{l+k 2 S, 2 ) 2 agrees perfectly 



with [TBI Unj . Importantly, it is characteristic for sound 
waves that localization is less pronounced at low energy. 
Indeed, for k£ —> at a fixed correlation ratio £ = cr/£ of 
order unity, the localization rate per wave length vanishes 
like [fc/i oc ] _1 ~ v 2 k£ — > 0. Consequently, the low-energy 
properties of the interacting quantum gas are not affected 
by localization. The underlying reason is that the disor- 
der is screened by interaction [TTJ — contrary to the case 
of noninteracting particles, where localization is stronger 
at lower energy |27j . 

In higher dimensions, localization is even less pro- 
nounced, with the localization length being exponentially 
large compared to the mean free path, if not infinite. 
Low-energy excitations are free to propagate over long 
times and large distances. The main effect of disorder 
then is to renormalize the excitation dispersion relation. 



DISORDER-MODIFIED DISPERSION 

The disorder-modified quasiparticle dispersion = 
efc+Aefe can be determined by applying standard Nambu- 
Green perturbation theory [55] to the relevant Hamilto- 
nian (12). Thanks to its perturbative structure H = 



jj(0) -|- H^ v \ it is straightforward to calculate the self- 
energy of the single-excitation Green function [26] . The 
self-energy's imaginary part provides the elastic scatter- 
ing rate, which vanishes at low energy just like the local- 
ization rate discussed above. From the self-energy's real 
part, we deduce the shift in the dispersion 



Ae fc 



v 2 a d 



d d <7 r ( v 



(17) 



with the kernel (P denotes the principal value) 



w 



(i) 

k(k+q) 



,(1) 12 
k(k+q)\ 

; + £fc+<7 



W 



(2) 

k(k+q) 



(18) 



Together with expressions (13) and (14), eq. (17) allows 



calculating the dispersion of Bogoliubov excitations in 
weak, but arbitrarily correlated disorder. 

In the hydrodynamic limit £ — > 0, where the healing 
length is shorter than both correlation length a and wave- 



length 2n/k, eq. (18) simplifies considerably, and eq. (17) 



reproduces eq. (28) of Ref.[14. In this regime, the exci- 
tation energy is reduced in all dimensions and for any 
value of ka. For low energies ka — > and smooth poten- 
tials a 3> £, the sound-velocity shift Ac/c = ~v 2 /(2d) is 
independent of the correlation details. 

We now change the point of view by taking the limit 
k — > 0, with an arbitrary correlation ratio £ = <r/£. In 
particular, this allows us to reach the case of (^-correlated 
disorder where a <C £, k~ x . The kernel ( 18 1 simplifies to 



ZOq 



q 2 e - (2 + q 2 e) COS 2 /3 



{2 + q 2 e) 3 



P = £{k,q). (19) 
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FIG. 1. (Color online) Sound-velocity correction, computed 
fro m (|17[ ) and (191, due to Gaussian correlated disorder, 
eq. ( |20p , with variance v 2 := V 2 /(gn) 2 < 1 as function of 
correlation ratio £ = cr/f. For most values, the sound veloc- 
ity is reduced. Dashed and dotted: universal limits for very 
smooth (£ 2> 1) and 5-correlated disorder ((" -C 1) respec- 
tively, as collected in Tab. [I] Inset: same data around the 
origin, showing the rapid departure from the leading-order 
estimate in d = 3 ISHTOl . 
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TABLE I. Universal limits of the speed-of-sound correction, 
computed from ( |17[ ) and ( |19[ ), for very smooth disorder (£ 3> 
1) and 5-correlated disorder <§C 1). See also footnote [] 



Fig.[T]shows the disorder correction ( 17 ) to the speed of 
sound, resulting from kernel ( 19 ) and a generic Gaussian 
pair correlation 



C d (ga) = (2^) d / 2 exp{- g V 2 /2} 



(20) 



The plotted curves can be expressed in closed form, but 
the details depend on the specific correlator and are not 
of general interest. In contrast, one finds universal be- 
havior for very small or very large £. The limit £ — > oo of 
a very smooth potential coincides, as it should, with the 
hydrodynamic limit Ac/c = —v 2 /(2d). In the opposite 
limit C — > of 5-correlated disorder, the correlator Cd(0) 
can be pulled out of the integral (17), which becomes el- 
ementary. The correction then scales as ( d , as shown in 
Fig. [I] with prefactors that are collected in Tab. |TJ [21] 
Notably, we corroborate the known result for (^-correlated 
disorder in d = 3 [8HI0] , the only case with a positive cor- 
rection. But our theory reveals that this estimate is of 
limited use because already a small correlation makes a 
large difference, as shown by the inset of Fig. [I] 



CONDENSATE DEPLETION 

Finally, we investigate the condensate depletion prop- 
erly speaking, namely the density of particles out of the 
mean-field condensate, Sn := L~ d J d d r^<5^i'(r)^*I'(r)). 
Inserting pj and rearranging terms, we find that the de- 
pletion density can be written 



Sn = r 

4nL d ^ 



k,k' 



akak'tik'-k 



nk'-k 



<7fc7fc' + 7fc'7k) 



nk'-k 
ak^k' 



{%1-k' +ltl- k ) - 2S ' 



fefe' 



(21) 



In principle, this expression is correct to all orders in 
V. One only requires the Fourier components ri\~ := 
[<E>(r) 2 ] fc = L~ d J2 q ®k-q® q of the deformed mean-field 
condensate density, as well as the Fourier components of 
its inverse, := n 2 [$(r) _2 ]fe. Up to order V 2 , they fol- 
low from the smoothing-theory results (|3| and Q. One 
also has to compute the Bogoliubov expectation values to 
the desired order. With the inhomogeneous Bogoliubov 



Hamiltonian (12 1 at hand, this is again a standard task 



in perturbation theory using the Nambu formalism |28j . 

Let us first check the homogeneous case 1^ = 0. Then, 
(21) shrinks to 6n<® = \ J (gi [a* ~ a^ 1 ] 2 , a well- 
known result [TS]. In d = 3, this evaluates to a deple- 
tion density <W 0) = [QV^tt 2 ^ 1 or equivalently to the 
relative depletion Sn^/n = Sfoal) 1 / 2 In d = 2, 
one finds Sn^ — [87r^ 2 ] -1 . The d = 1 integral is in- 
frared divergent, consistent with the fact that zero-point 
fluctuations prevent homogeneous ID BECs. Cutting 
off the integral at some value a = £ki& <C 1, with 
kin of the order of the inverse system size, one finds 
6nW = (2 In 2 - 2 - lna)/(2 v / 2<), up to order a. 

Now we evaluate the disorder-induced depletion by ex- 
panding all contributions to ( pTj ) to second order in V. 
Upon taking the ensemble average, terms of order V av- 
erage to zero, and Sn = Sn^ +Sn^ +0(v 3 ). Each of the 
second-order terms is proportional to the correlator (151, 



and we can collect all contributions into a single kernel: 



fe( 2 ) = v 2 Sn^ 



d d q 
(27r)° 



a d C d (qa)G d (qO- 



(22) 



Details of this derivation will be given elsewhere (30) . 
Here we note that the relative depletion in units of u 2 , 
A(C) = Sn^ /v 2 6n^ is only function of C = cr/f. [31] 
This quantity is plotted in Fig. [2] As for the sound ve- 
locity correction, details depend on the specific correla- 
tion. In the (5-correlated limit C - > 0, Fig. [2] shows the 
generic scaling A(£) = (3d( d Cd(0) with numerical coef- 
ficients j3i ps 0.235 (weakly dependent on the cutoff a), 
/3 2 ~ 0.135, & ps 0.160. In the limit C ^ oo of a very 
smooth potential, we find A — > Gd(0) with C?3(0) = 3/8, 
G 2 (0) = 0, and Gi(0) = -1/8. In all cases but the last, 
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FIG. 2. Disorder-induced condensate depletion (22 1, rela- 
tive to the clean value and in units of disorder strength v 2 , 
as function of the correlation ratio f = cr/£ for the Gaus- 
sian correlation poj. Dashed and dotted: universal limits as 
collected in Tab. IHI 



in all dimensions. Moreover, we calculate the disorder- 
induced quantum depletion, which proves to be much 
smaller than the previously known mean-field conden- 
sate deformation. We conclude that our theory should 
fare very well in describing the excitations of disordered 
interacting bosons, especially in dilute cold gases, where 
the study of well-controlled disorder in earnest has just 
begun [liHZg. 
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TABLE II. Universal limits of the disorder-induced conden- 
sate depletion, relative to the clean value and in units of v 2 
as plotted in Fig. [2] for very smooth disorder (£ 3> 1) and 
5-correlated disorder <C 1). For the latter case, the numer- 
ical coefficients are Pi ~ 0.235 (for a — 0.01), /?2 ~ 0.135, 
/3 3 w 0.160. 



the condensate depletion due to disorder is positive, as 
expected. As shown in Fig. [2j also in d = 1 the depletion 
is positive for most values of £. The curve only crosses 
over to negative values for such a large value C = °7£ 
depending on the cutoff a, that the correlation length a 
has to be comparable to the system size, which is not the 
regime of present interest. 

In all cases, the combined depletion due to interaction 
and disorder reads Sn = <5n(°)[l + w 2 A(C)], with A(C) at 
most of the order of unity. Clearly, the fractional deple- 
tion induced by the disorder, 5rS 2 ^ /n = (dnA°) / 'n)v 2 A(£) 
is a factor 8n^ /nCl smaller than the mean-field con- 
densate deformation, which is of order v 2 . In hindsight, 
this result is rather plausible: the primary effect of the 
external disorder potential is merely to deform the con- 
densate. The depletion of the condensate itself is a sec- 
ondary scattering effect, mediated by the weak repulsive 
boson interaction, and therefore considerably weaker. 

In conclusion, we report substantial progress in the 
analytical description of interacting, condensed bosons 
in correlated disorder of any dimensionality. We de- 
rive the fundamental Bogoliubov Hamiltonian for exci- 
tations. This determines a wealth of (thermo-)dynamic 
quantities, out of which we calculate the sound velocity 
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